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Abstract 



Lattice work, exploring the Higgs mass triviality bound, seems to indicate that a strongly interacting 
scalar sector in the minimal standard model cannot exist while low energy QCD phenomenology seems to 
indicate that it could. We attack this puzzle using the 1/N expansion and discover a simple criterion for 
selecting a lattice action that is more likely to produce a heavy Higgs particle. Our large N calculation 
suggests that the Higgs mass bound might be around 850 GeV, which is about 30% higher than previously 
obtained. 
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The recent quantitative results that emerged from the investigation of triviality in the A((/> 2 ) 2 theory 
on the lattice [1] have given an upper bound on the ratio of the Higgs mass, tuh, to the weak scale, 
/tt = 246 GeV, of about 2.5-2.7. In the hypothetical world of QCD with just the up-down doublet of 
quarks with zero masses, which is described by the same effective theory at low energies, one believes, on 
phenomenological grounds, that the above ratio would be significantly larger, equal to about 6-8. This is 
admittedly an extreme situation because the sigma enhancement in n n scattering doesn't quite qualify to 
be called a resonance, and, it comes accompanied by a host of other resonances, representing all kinds of 
nontrivial "cutoff effects". Still, there is a large gap between 2.7 and 6 and we are left with essentially two 
alternatives: either the QCD example provides a totally unreasonable estimation for the bound, or, we admit 
that the lattice result, while solidly established within a subclass of actions, is relatively far from a realistic 
estimate for the bound. We would have much less cause for worry if the lattice bound for the ratio came out 
to be about 4 (this would correspond to to# ~ 1 TeV) or more [2]. 

This issue could be addressed in principle by an extensive search of all possible lattice actions. In 
practice one needs some theoretical indication as to what kind of regularized action is likely to produce a 
higher bound. If the physical picture behind the particular form of the suggested action is reasonable, an 
investigation of its consequences would help resolve the above puzzle by either refuting the existing lattice 
bound or by strengthening the contention that the QCD example is grossly exaggerated. 

The purpose of the present note is to make a proposal for such a generalized action and sketch the 
theoretical basis for it. The basic point of the proposal is that one should put on the lattice nonlinear 
(9(4) models that have a naive continuum limit with freely adjustable four derivative couplings and study 
them nonpcrturbatively. We solve the model in the large N limit of O(N) and find the region in coupling 
space where the highest bound is expected; we test the extent to which the results are dependent on the 
rcgularization and extract the basic features that are generic. The results can be explained by a rough 
physical picture as follows: When the regularized model is nonlinear one has to think about the Higgs 
resonance as a loosely bound state of two pions in an I = 0, J = state. Pions in such a state attract 
because superposing the field configurations corresponding to individual pions makes the state look more like 
the vacuum and hence lowers the energy. The four derivative terms in the action can add to this attraction; 
if the attraction is increased to a certain level the bound state will become massless and stable. This occurs 
while the system is still relatively deep in the broken phase; hence the weak scale f n is nonvanishing. The 
appearance of this massless bound state signals a new kind of phase transition; the ordinary phase transition 
occurs when both f n and mj? vanish and leads to a scale invariant theory at the transition - here we still 
have a scale and therefore scale invariancc seems to be spontaneously broken, with the massless boundstate, 
the Higgs particle, playing the role of the expected dilaton. This phenomenon is reminiscent of the results 
obtained in three dimensions by Bardeen, Moshe and Bander [3]. As far as the Higgs mass bound goes, it 
is clear that we have to try to be as far away as possible from this new transition point; we wish to induce 
enough repulsion between the pions to delay the creation of the bound state as much as possible. Of course, 
the interaction is still attractive overall, because only the order p 4 term in the interaction can be tampered 
with; the order p 2 and p log(p ) terms are fixed by current algebra and are attractive [4]. Moreover, in 
the critical regime the model still is governed by the A((/> 2 ) 2 trivial fixed point and is expected therefore to 
generate a one particle state corresponding to the a field component. While it is quite possible that the 
above described scale invariance breaking transition cannot actually materialize at N = 4 (we shall see that 
at N = oo it does appear in all regularizations that we tested), the basic trend for the magnitude of the 
triviality bound must be the same. 

Assuming global 0(4) symmetry and that the leading cutoff effects are represented by dimension six 
operators, one arrives, after eliminating couplings that won't lead to observable effects in the S-matrix to 
leading order in the inverse cutoff, to a space of actions that depends on 4 parameters [5]. One may try to 
induce the dependence on the associated renormalizcd operators by simply writing identical forms for terms 
in the bare action. However, we are looking for a bound on the physical four point coupling and it is likely 
that it obtains in the limit that the bare four point coupling is taken to infinity. When the bare coupling 
becomes infinite the bare model is reduced to a nonlinear model; this has the effect that the naive forms of 
the dimension six operators in terms of bare linear fields become trivial. The parameterization in terms of 
linear fields and operators of naive dimension six breaks down completely at the point where the bare theory 
becomes nonlinear (we have tested this explicitly in the linear model, within 1/N, and with a particular form 
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of a bare dimension-six term in the Lagrangian). So we were expecting to be left with three freely adjustable 
parameters (within some range) but we ended up with only one. This is the case that has been investigated 
nonperturbatively until now [1] . To include the effects of dimension six operators at physical scales we need 
to put in operators of higher order in derivatives in the bare, nonlinear, theory that describes the physics at 
cutoff scales. There are exactly two additional terms that one can write down if one restricts the number of 
derivatives to four [4]; thus the right number of parameters is obtained. This time, unlike in the linear case, 
the effect of the bare couplings has no reason to disappear. 

It is well known that the 1/JV expansion, at least to leading order, doesn't have to work well quantita- 
tively simply because there is a factor of N + 8 in the coefficients of the /? function and therefore N = 4 is 
badly approximated by N = oo [6-8]. We shall say more about this later; for the moment keep in mind that 
we are interested mainly in relative effects, that is, in how the bound varies with the action, and less in the 
absolute magnitude of the bound. To get the latter one would have to resort to numerical means. 

If a strongly interacting scalar sector could be produced on the lattice its investigation would require 
some nontrivial extension of the tools that were used until now. However, if we only wish to establish the 
existence of a strongly interacting regime, we can do that indirectly, while restricting our attention just to 
weak couplings. For example, a strongly coupled Higgs sector was excluded on lattices with a pure nearest- 
neighbor coupling, by showing that, for a ratio mn/f-K of only 2.5, one already had to^/A equal to 0.5, 
leaving no "room" below the cutoff A for further growth in the coupling; if we could show that a different 
action would, at the same value of ma/fir, give W///A equal to 1/10 say, our case would be made. Therefore, 
our subsequent study is restricted to the weakly coupled sector. Generically, we expect there that 

^«C(6)cxp[-16^ 2 / 2 /(A^4)L (1) 

where C(b) is a function of the additional parameters we introduced in the action. We aim to show that 
by varying b C(b) can be significantly lowered below the value it had when no four derivative terms were 
included (this is the only case that has been investigated to date). Our finding is best summarized in a 
qualitative self-explanatory diagram that we checked for a class of Pauli-Villars regularizations and for the 
F4 lattice regularization (the generalization of the calculations to the hypercubic lattice should be a simple 
exercise) . 

We first look at the models with Pauli-Villars, or higher derivative, regularization. We choose to replace 
the ordinary propagator by l/[p 2 (l + (p 2 /M 2 ) n )] where M is a "bare" cutoff (as opposed to A the physical 
cutoff to be defined below) and n > 3 is an integer. The limit n = 00 corresponds to a sharp momentum 
cutoff and is expected to be singular. For any finite n the bare propagator has n "ghost" poles all situated 
on the unit circle in the complex s/M 2 plane. We shall define the physical cutoff, A, as the square root of 
the modulus of the location of the closest complex pole in s in the full pion propagator (at leading order in 
1/N). The action, written in Euclidean space, is: 



S = 



(2) 



where g(-d 2 ) = (-<9 2 )[1 + (-<9 2 /M 2 )"], and the field is subjected to the constraint (jr = N/3. The partition 
function is Z = f[dcp] exp[— S]. 

Using the technique employed in ref. [9] we find the phase diagram of the model paying particular 
attention to the competition between various saddles. We set 

$= V~Nv + vH(x) + ir(x) (3) 

with v — v/v, v — \v\, v-tt = 0, / H(x) = 0, f ir(x) = 0. We define the pion wave function renormalization 
constant Z v for the field n as usual. The pion decay constant f n can be easily gotten from v and Z v . 
It turns out that it is more convenient, for the broken phase, to parameterize the theory by Z^: 

_ ( 4tt\ 2 nsin(27r/n) Z„ - 1 
1 \M 2 ) 2tt zl (n ~ 1)/n 



3 



3 ion Wave Function Renormalization Constant 




Fig. 1 : Summary Diagram 

Qualitative Features Common to all Models Investigated 



where < Z n < 2(n — l)/{n — 2). In this range, b\ increases monotonically with and at the larger 
end-point we have dbi/dZ n = 0. The portion of the phase diagram we are interested in does not depend 
on 6 2 because giving an expectation value to the O(N) singlets associated with 6 2 would break Lorentz 
invariance. The other coupling is best replaced by the unitless weak scale, F = / W /A. The physical cutoff, 
A, is given by A = M/zl^ 2n \ From the large N saddle point equations we obtain the expectation value v 
which is related to f„ by v 2 — Z^f 2 resulting in: 



N ~\m) P * 



nsin(7r/n) 



(5) 



There is a second order transition line at F — and, on it, at the end of the allowed range of Z„, there 
is a tricritical point where the second order line turns into first order. At the tricritical point one obtains a 
massless pole in the two point correlation function of the trace of the energy-momentum tensor. We shall not 
be concerned with the region beyond the tricritical point. The calculation of ran proceeds in the standard 
manner; it turns out that one has to find the complex zeros of the determinant of a four by four matrix. 
Since we restrict our attention to the perturbative regime, i.e. F, p << 1 we can replace the quantity ran 
by the technically more accessible real quantity ran, defined from the smallest positive zero in s = — p 2 of 
the real part of the determinant by ran — s/s. We are really interested in Mr — ran/ A. If we also take 
Mr « 1 of the same order (up to logs) as F, the fourth row and column of the matrix decouple and the 
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dependence on 6 2 disappears. We need the expression for the "bubble" diagram 
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\og{p 2 /M 2 ) 



(6) 



+ 0(p 2 /M 2 ) 



which holds for any n < 00 (for n = 00 the correction term goes as \Jp 2 /M 2 and represents a non-local term 
in the effective action as expected when employing a sharp momentum cutoff). One has then: 



M 2 R w exp 



_2_ UirFY 
'N \~Mr~) 



C(Z 7T ) = exp 



n — 1 
2n 



1 - 



1 - Z v 



ntan(7r/n) 1 - Z T (n - 2)/(2(n - 1)) 



(7) 



There is no four derivative term when Z„ = 1; the minimal C obtains at Z„ = which is as far away from 
the tricritical point as possible: 



C(Z„ = 1) 
C{Z V = 0) 



exp 



1 



2n J ntan(7r/n) 



< exp(0.5) w 1.65. 
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We have obtained some increase in the cutoff but not a factor of about 5 as we would have liked. 
However, we know from the work of Einhorn [7] that even at bi = this particular large N model will have 
a strongly interacting sector; therefore it may be that there is no room left to dramatically increase the 
cutoff even further, ft is also clear that if it turned out that the lattice actions that were so far investigated 
really correspond to an action with b\ > in Pauli-Villars regularization, a much larger variation could be 
obtained because C diverges at some finite, positive b\. b\ > means that the extra term has induced more 
attraction in the I,J = channel and the divergence in C occurs at the tricritical point where the resonant 
bound state becomes massless and stable. 

We turn now to the lattice. Since our objective here is to provide a useful indication for what kind of 
actions would be instructive to study by nonperturbative means, we pick the lattice regularization that is 
best suited for this purpose, namely we put the theory on an F4 lattice. This lattice has a larger symmetry 
group than the more commonly used hypercubic lattice which ensures that no Lorentz invariance breaking 
four derivative terms appear in the (naive) continuum limit of its lattice actions. The F4 lattice can be 
viewed as a hypercubic lattice (Z 4 ) with unit lattice spacing from which all the sites with an odd sum of 
coordinates have been erased [5]. We shall denote sites on F4 by x,x',x", and links by < x,x' >,l,l'. The 
action is chosen as 

S = - 2N/3 ®( x ) ■ ®( x ') - N @i zZ [®( x ) ' ®( x ') 



48 Z-~i 



48 



E 



$(x) • <S(x') 



(9) 



Here the field is constrained by $ 2 (x) = 1. To express the action in a form closer to the one we employed 

before we rescale the field <j) — ^/6-/V(/3 + fh + #2)$ ( we om y consider the region /3 + Pi + P2 > 0). The 
action has the same form as above (eq. (2)), up to terms that contain fourth order derivatives, if one 
identifies: 



1 Pi 



■02 



48 (p + Pi+P 2 ) 2 '- 



b 2 



12 

01 



(10) 
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Again 62 plays no role and one may set (3± to 0. The situation is essentially the same as in the Pauli- 
Villars case. b\ can be traded for the pion wave function renormalization constant, 

where Z v varies between and 2. The standard non-linear action corresponds again to Z w = 1 and, on the 
line where the pion decay constant F vanishes, = 2 corresponds to the tricritical point. F is given by: 

-rf = = r , r = / ——= 0.13823. (12) 

N Z n J k g{k) 

The momentum integral is over the F 4 Brillouin zone and a factor of 1/ (2ir) A has been absorbed in the 
definition of J, ; the numerical value is taken from rcf. [5]. 

We proceed now to find Mr defined exactly as before; for this we need the lattice "bubble" which can 
be easily evaluated using the results of ref. [5] (but note a missing minus sign in the value of n in equation 
(4.12) there): 

B F , (p) = - lo §(P 2 ) + c i + C 2P 2 + (P A lo §P 2 )' 

(13) 

ci = 0.0466316, c 2 = 0.0005497. 

Unlike in the Pauli-Villars case, the cutoff (defined for example by a certain amount of violation of Euclidean 
rotational invariance in the full pion propagator) is independent of the couplings in the action. 

This time we shall be interested to go beyond the 0(1) term and that's why we evaluated the next 
correction in p to the bubble (13). Ignoring this correction for the moment we get the usual formula for Mr 
(eq. (7)), but the coefficient C is now given by, 



C(Z 7T ) = exp 

leading to: 



8n ( Cl - r °T3zj2 



(14) 



g^=«p(*Ag)« 4.521. (15) 

We see that we have obtained a much larger variation. 

In practice, for vectorization purposes, it would be better to have an action that involves only nearest 
neighbors on the F4 lattice. Nothing in what we have said until now changes if we replace the /?2 term in 
the action by 

n^y, E [(*(*) • • • ( 16 ) 

x <W> 

1,1'ni^l Inx'jtHI, I'n.'VH, x,x',x" all n.n. 

To check sensitivity to higher order terms in F 2 and/or M R , while still neglecting & 2) we go to the next 
higher order in F, p, Mr when we find the Higgs mass. It turns out that the numerical correction is small (of 
the order of 1-2 % ) even for Mr ss 1. For the usually studied action with b\ and 62 = we can check how well 
the 1/N expansion works, on the quantitative level, in the region of main interest, namely ttir/ f n w 2.0 — 2.5; 
in this region we have numerical data at -/V = 4 (see last reference in [1]). We find agreement to within 
25 percent for the Mr/F ratio as a function of Mr] the large N results are systematically larger than the 
N = A ones but the relative change in the ratio Mr/F when Mr is varied is more accurately predicted. 
This is enough to convince us that our large N results ought to be taken seriously as an indication for what 
is going to happen in the real system. For example, we would predict that if we limit cutoff effects by the 
reasonable [5] requirement Mr < V2/2 we find, with N = 4: 

Mr 2ir 

~ET < I (17) 
b y/log [V2C(Z V )] 
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This gives for the nearest neighbor case (Z^ = 1) Mr/F =3.1 while for the extreme case (Z v = 0) one finds 
Mr/F = 4.0. For N = 4 the numerical result is smaller by about 20%. We therefore are led to conjecture 
that for the simplest case with j3\ — the bound would increase by 30% to about 850 GeV . 

We should comment here on the possible relationship between this note and the work based on effective 
Lagrangians that was carried out in the context of the investigation of the scalar sector of the standard model 
by several groups recently [10]. Although in both cases one employs a nonlinear Lagrangian of identical forms 
(at least for N = 4) there is a fundamental difference that ought to be stressed: In the effective Lagrangian 
approach one uses renormalized perturbation theory as defined by Weinberg [4] to summarize all relations 
between soft pion amplitudes to any finite order in the external momenta of the pions; one then tries to 
extrapolate [10] these results by several schemes, all attempting to approximately impose unitarity and 
crossing, to the region of energies dominated by the tails of the lowest resonances. At no point does this 
approach consider the question of whether the physics so defined can be consistently extended to the cutoff 
(see second reference in [4] for some discussion). Technically speaking, chiral renormalized perturbation 
theory, very much like ordinary perturbation theory, does not provide the relation between the parameters 
in the bare Lagrangian (i.e. cutoff scale physics) and the parameters that describe low energy properties; 
it only provides a very large number of constraints between all the low energy physical processes. Thus, 
questions having to do with triviality are completely outside the reach of these approaches. In more intuitive 
terms, we are taking here very seriously a model defined by a chiral Lagrangian, all the way to the highest 
energies, and seek for constraints on the allowable low energy effects that result from the requirement that 
they occur below the cutoff. If the cutoff is very high relative to the lowest resonance we are bound, by the 
general principles of the renormalization group, to view this as just another regularization of A(</> 2 ) 2 and 
hence a perfectly reasonable framework for investigating consequences of triviality. 

We have not discussed in any detail the term with coupling 6 2 - We have seen that it has no effect in 
the perturbative regime; however, when terms beyond the equivalent of the 0(p 2 ) correction to the bubble 
are taken into account, 62 comes into play and can become important. The structure of the term indicates 
that 62 is the right dial to tune if we wish to make our theory look more like QCD. This suggests that it 
might help to make the Higgs mass bound even higher. But in the present note we were interested mainly 
in estimating the Higgs-mass triviality bound when nothing special (technicolor-like) has been imposed on 
the theory. Further discussion of the 62 term would take us too far afield. We shall also omit any discussion 
of subleading terms in the 1/N expansion; particularly interesting would be the corrections to Z w and to 
the phase structure in the vicinity of the tricritical point. These topics and a comprehensive analysis of 
Renormalization Group flows in the models under consideration will hopefully be addressed in future work. 

We should also comment here on the fact that it is well known that once the scalar sector is strongly 
interacting, looking at rriR rather than mn, is a bad approximation due to the large width [11]. This effect 
leads to a saturation of the growth of mn much before tur stops growing and can be best seen in the 1 /N 
expansion, as shown by Einhorn [7]. It is important to realize that this sort of bound has nothing to do 
with triviality; it appears even when the issue of triviality is totally ignored, for example, when the bare four 
point coupling in the linear model is allowed to become negative and the theory has a tachyonic instability. 
Our calculations can all be done for a relatively narrow Higgs, and all they aim to show is that a strongly 
interacting scalar sector in the minimal standard model has not yet been convincingly ruled out. The relation 
between vhr and mn alluded to above may very well hold in such a sector, if only we could show that one 
exists. 

In summary, we propose that the scalar sector of the minimal Standard Model can be made to have a 
stronger renormalized self coupling by introducing, at the cutoff level, four-derivative couplings that induce 
some repulsion between Goldstone Bosons. Within the 1/N expansion this increases the Higgs mass bound, 
compared to the usual lattice result, by up to about 30%. Obviously it would be extremely interesting to 
study, by numerical means, whether this conclusion persists for the physically relevant case of N = 4. This 
would help resolve the puzzle presented at the beginning of this Letter. We hope to address this question in 
the future. 
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